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Abstract 

We solve for the complete space of linearized deformations of the Klebanov-Strassler back- 
ground consistent with the symmetries preserved by a stack of anti-D3 branes smeared on 
the S s of the deformed conifold. We find that the only solution whose UV physics is consis- 
tent with that of a perturbation produced by anti-D3 branes must have a singularity in the 
infrared, coming from NS and RR three-form field strengths whose energy density diverges. 
If this singularity is admissible, our solution describes the backreaction of the anti-D3 branes, 
and is thus likely to be dual to the conjectured metastable vacuum in the Klebanov-Strassler 
field theory. If this singularity is not admissible, then our analysis strongly suggests that 
anti-D3 branes do not give rise to metastable Klebanov-Strassler vacua, which would have 
dramatic consequences for some string theory constructions of de Sitter space. Key to this 
result is a simple, universal form for the force on a probe D3-brane in our ansatz. 



1 



Contents 



1 Introduction 

1.1 Strategy and Results 

2 Perturbations around a supersymmetric solution 

2.1 The First Order Formalism 

2.2 Papadopoulos-Tseytlin ansatz for the perturbation 

2.3 The zero-th order solution: Klebanov-Strassler .... 

2.4 £j Equations 

2.5 <p % Equations 

2.6 Imaginary (anti)-Self Duality of The Three Form Flux 

3 The Force on a Probe D3 Brane 

4 The Space of Solutions 

4.1 Solving the £j Equations for X\ — 

4.2 Solution for X\ ^ 0, UV and IR expansions 

4.3 Solving the (f) % Equations 

4.3.1 The Space of Solutions 

4.3.2 IR behavior 

4.3.3 UV behavior 

4.4 The Zero Energy Condition and Gauge Freedom . . . 

5 Boundary Conditions for Three-Branes 

5.1 BPS Three-Branes 

5.2 Anti-BPS 3-branes in Anti-Klebanov-Strassler .... 

5.3 Anti-BPS 3-branes in Klebanov-Strassler 

6 Constructing the Anti-D3 Brane Solution 

6.1 UV IR, X\ is non zero 

6.1.1 IR Singularities 

6.2 IR ^ UV 

7 Discussion 
A Conventions 

B Behavior of £ 

B.l IR behavior of £ 

B.2 UV behavior of £ 



1 Introduction 

Over the past few years the analysis of metastable vacua of supersymmetric field theories has begun 
to play an important role in phenomenological model building, as these vacua can bypass many 
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of the familiar problems associated with old-fashioned dynamical supersymmetry breaking (DSB). 
Such meta-stable vacua were first found for massive M = 1 SU(Nc) SQCD in the free-magnetic 
phase [lj , and soon thereafter in quite a large number of similar theories. 

An obvious question posed by the ubiquitousness of metastable vacua in SQCD theories is 
whether such vacua also exist in the string-theory configurations whose low-energy physics yields 
these theories. These constructions involve D-branes and NS5 branes in type IIA string theory, 
and are generically referred to as MQCD [2J. Another obvious question is whether such metastable 
vacua exist in the type IIB supergravity backgrounds such as [3] that are dual to some of the SQCD 
theories via the gauge-gravity (AdS-CFT) correspondence. 

For MQCD theories one can construct putative non-supersymmetric brane configurations whose 
low-energy physics matches that of the SQCD metastable vacuum. At g s = these brane config- 
urations have the same asymptotics as the configuration describing the MQCD supersymmetric 
vacua. However, for g s ^ these non-supersymmetric configurations do not have the same asymp- 
totics as the MQCD supersymmetric branes, and hence do not describe a metastable vacuum 
of a supersymmetric theory (where supersymmetry is broken dynamically), but rather a non- 
supersymmetric vacuum of a nonsupersymmetric theory [4]. The fact that the putative MQCD 
non-supersymmetric and supersymmetric configurations are not vacua of the same theory can also 
be ascertained by computing the tunneling probability between them: for g s ^ this probability 
is zero because the branes differ by an infinite amount, while for g s = this probability is zero 
because the branes have infinite mass. 

The key reason for which the susy and non-susy brane configurations do not have the same 
asymptotics comes from the phenomenon of brane bending: essentially all MQCD brane descrip- 
tions of four- dimensional theories contain D4 branes ending on NS5 branes; since the end of D4 
brane is codimension 2 inside the NS5 brane, it sources an NS5 worldvolume field that is loga- 
rithmic at infinity, and that causes the NS5 brane to bend logarithmically. For the MQCD dual 
of M = 2 SQCD, this bending encodes the logarithmic running of the coupling constant with the 
energy, characteristic of four-dimensional gauge theories J2J. The D4 branes in the deep infrared 
of the susy and non-susy MQCD brane configurations end on the NS5 branes in different fashions, 
and source fields of different kinds. Since these fields do not decay at infinity, but grow logarith- 
mically, the tiny infrared difference is transposed into a log-growing difference in the UV, which is 
responsible for the susy an non-susy configurations not being vacua of the same theory. 

It is also interesting to note that for three-dimensional gauge theories engineered using D3 
branes and NS5 branes in type IIB string theory, this phenomenon does not happen. The D3 
branes end in codimension-three sources on the NS5 branes, which source fields that decay like 1/r 
at infinity. This bending encodes the running of coupling constants in three dimensions, which is 
linear in the inverse of the energy. Since the difference between the fields of the non-supersymmetric 
and supersymmetric brane configurations decays at infinity, the two do correspond to vacua of the 
same theory, and the non-supersymmetric configuration is really a 3D MQCD metastable vacuum. 

Thus, the mechanism that ensures the absence of metastable vacua in four- dimensional MQCD 
appears to rely on just one ingredient: the existence of logarithmic modes that are sourced dif- 
ferently in the infrared for the BPS and non-BPS configurations. As these modes are present in 
four- dimensional but not in three dimensional gauge theories, it is only four- dimensional MQCD 
which does not have metastable vacua, three-dimensional MQCD does have such vacua. 

In this paper we will analyze whether a similar mechanism is at work in gravity duals of M = 1 
four- dimensional gauge theories. As the coupling constants of these theories run logarithmically 
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with the energy, all their gravity duals must contain fields that depend logarithmically on the 
AdS radius, at least asymptoticallyS. It is possible that a putative metastable configuration in 
the infrared couples to these logarithmic modes in different manner than the supersymmetric 
configuration; hence the non-BPS solution would differ in the ultraviolet from the BPS solution 
by one or more non-normalizable modes. If this indeed happens, then the non-BPS solution will 
not be dual to a metastable DSB vacuum of the supersymmetric field theory. Instead, it will 
describe a non-supersymmetric vacuum of a non-supersymmetric theory obtained by perturbing 
the Lagrangian of the supersymmetric theory. As a consequence, supersymmetry breaking will not 
be dynamical but explicit. 

The best-known candidate for the IIB gravity dual of a metastable gauge theory vacuum has 
been proposed by Kachru, Pearson and Verlinde (KPV) [5], who argued that a collection of P 
anti-D3 branes placed at the bottom of the Klebanov-Strassler (KS) solution with M units of 
3-form flux can decay into a BPS solution corresponding to the KS solution with M — P D3 
branes at the bottom of the throat. In particular, if P — M this collection of anti-D3 branes 
could decay into the smooth, source-free KS solution. The argument presented by [5] was in 
the probe approximation and it important to extend this to a backreacted supergravity solution. 
Finding the backreacted solution would allow one to establish whether this configuration is indeed 
a metastable vacuum of a supersymmetric theory, or it is rather a non-supersymmetric vacuum of 
a non-supersymmetric theory. This necessity is particularly evident, given that in MQCD there 
exist non-backreacted brane configurations that are candidates for meta-stable vacua [6j [7J 0] but 
in fact break supersymmetry explicitly [I] once backreaction is taken into account0. 

Note that this question is in no way settled by the results currently in the literature. For 
example the ultraviolet perturbation expansion around the Klebanov-Tseytlin [UJ solution pre- 
sented in [12] assumes a-priori that the anti-D3 branes source only a normalizable mode, and 
uses the consistency of the result (in particular they use the fact that the force on a probe D3 
brane is not inconsistent with that calculated in [13]) to argue that such an assumption is correct. 
However as we will show in Section [31 all SU (2) x 577(2) x Z 2 -invariant first-order perturbations 
around the Klebanov-Strassler solution, whether normalizable or non-normalizable, give a force 
on a probe D3 brane that is either zero or has the precise behavior predicted in [13]. Hence this 
force calculation alone will not determine the normalizability or non-normalizabilty of the modes 
sourced by an anti-D3 brane. In addition, we will point out that subleading perturbations of the 
(singular) Klebanov-Tseytlin geometry do not in fact correspond to perturbations of the (regu- 
lar) Klebanov-Strassler geometry and as such the Klebanov-Tseytlin geometry is an inappropriate 
arena for studying these issues. 

One can argue that because an anti-D3 brane is a point source in a six- dimensional space, the 
fields sourced by such a brane should decay at infinity as 1/r 4 , and thus could not correspond to 
non-normalizable modes. As we will show in Section I6\T1 by explicit calculations, this intuition is 
not correct even under the weakest of assumptions about the backreacted solution sourced by anti- 
D3 branes, which must decay at least as strongly as 1/r 3 . Furthermore, as we have argued above, 
both the Klebanov-Strassler solution, as well as other backgrounds dual to four-dimensional gauge 

the well-known Klebanov-Strassler solution [3] such a field is given by the integral of the NS B-field on the 
2-cycle of the conifold. 

2 Worldsheet evidence in support of this MQCD analysis was given in [8] and furthermore, both the MQCD dual 
of the configuration of [5] (analyzed in [9]) as well as other duals of KPV-like configurations 10 , suffer from the 
same problems in the UV, at least in certain regime of parameters. 
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theories whose coupling constants run logarithmically with the energy, contain supergravity modes 
that behave asymptotically as logr. If the anti-D3 brane in the deep infrared couples to such a 
mode, then the backreacted solution would differ from the BPS solution by a non-normalizable 
mode, and hence would not correspond to a metastable DSB vacuum. 

It is also interesting to note that the only known fully-backreacted supergravity solution cor- 
responding to a DSB metastable vacuum, constructed by Maldacena and Nastase [H] . describes 
a vacuum of a 2+1 dimensional gauge theory, for which none of the arguments presented above 
applies. As we have discussed above, these theories also have MQCD metastable DSB vacua, and 
their coupling constants are linear in the inverse of the energy, and hence do not correspond to 
bulk fields that grow in the UV and that could amplify the tiny infrared difference between the 
BPS and non-BPS solutions into a non-normalizable UV mode. 

Our purpose in this paper is to establish the nature of the fields that are sourced by anti-D3 
branes placed at the bottom of the Klebanov-Strassler solution, and to see whether these branes 
give rise to normalizable or non-normalizable modes. To do this in generality one would need to find 
the fully-backreacted solution corresponding to anti-D3 branes in the warped deformed conifold. 
Since these branes are point sources in the internal space, and since moreover they polarize (by 
the Myers effect [TJ]) into D5 or NS5 branes [5J, the fields of such a non-supersymmetric solution 
would depend on three or four variables. Finding such a solution is clearly beyond hope with 
current technology. 

One way to simplify the problem is to smear the anti-D3 branes on the £ 3 of the deformed 
conifold, so that the resulting backreacted solution preserves the £77(2) x £77(2) x Z2 symmetry of 
the original solution. The ansatz for the most general solution preserving this isometry was written 
down by Papadopoulous and Tseytlin [IB] , who reduced the supergravity equations of motion to 
second-order ordinary differential equations for eight scalar functions of the "radial" variable r of 
the KS solution. 

One might object that this smearing washes down important features of the anti-D3 branes that 
give the putative metastable vacuum. Indeed, as discussed before, when these anti-D3 branes are 
together they acquire extra dipole moments by polarization, and hence couple to more supergravity 
fields than the smeared D3 branes. Furthermore, because of the broken supersymmetry one can 
argue that the anti-D3 branes cannot remain smeared on the £ 3 , but will eventually condense into 
a single-center configuration. This condensation however takes place over a time that can be tuned 
to be parametrically large, and hence does not affect the fact the that smeared anti-D3 branes give 
rise to consistent supergravity solution. Furthermore, as the smeared anti-D3 branes couple to less 
fields than the localized polarized ones, their potential for sourcing a potentially-dangerous mode 
that is non-normalizable in the ultraviolet will be much lower than for the localized and polarized 
anti-D3 branes. Hence, if the smeared branes give rise to non-normalizable modes, the localized 
branes will very likely do the same; however, if the smeared branes do not source non-normalizable 
modes, this does not exclude at all the possibility that the localized polarized branes will source 
them. 

Even after smearing the anti-D3 branes, solving the system of eight coupled nonlinear differen- 
tial equations is no easy task. To make the general problem tractable one must use the fact that 
the fields sourced by the anti-D3 branes are subdominant compared to the fields of the background, 
at least in the region far away from the tip. Hence, one can study these fields by performing a 
first-order perturbation expansion around the supersymmetric KS solution. 

Fortunately, as shown in detail in [IT] , the second-order equations satisfied by the perturbations 
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of the eight scalar functions of the PT ansatz around the KS solution, factorize into sixteen first- 
order equations of which eight form a closed system. Note that this factorization does not throw 
away any of the modes: both the eight second-order equations and the sixteen first-order ones 
have sixteen integration constants. In fact, as we will show in Section l4~4l one of these integration 
constants is a gauge artifact and another one must be fixed to zero in order for the solutions to 
the sixteen equations to give a solution to Einstein's equations. Hence, the class of perturbative 
solutions that we construct depends on fourteen physically-relevant integration constants! 

Before beginning it is worth discussing the implication of the potential existence of a non- 
normalizable mode sourced by anti-D3 branes at the bottom of the Klebanov-Strassler solution. 
Such anti-branes are part of the "staple diet" of string phenomenology and string cosmology 
constructions: they are a crucial ingredient in the KKLT construction of de Sitter space [H] and 
in the KKLMMT model for string theory inflation [13]. These constructions involve placing anti- 
D3 branes in a Klebanov-Strassler throat that is glued to an ambient Calabi-Yau. The anti-branes 
lift the energy of the AdS vacuum giving rise to de-Sitter vacua and can create a potential for a 
probe D3 brane that drives inflation. Furthermore, one generally assumes that the extra energy the 
branes bring is parametrically redshifted in the KS throat, and can give rise to mass hierarchies. 

A non-normalizable mode would "climb" up all the way to the end of the Klebanov-Strassler 
throat, and will probably affect the gluing of this throat to the ambient Calabi-Yau. However, as 
this gluing is not understood, it is premature to say whether it will be negatively affected by the 
non-normalizable mode. Nevertheless, if a non-normalizable mode is present, its energy will not 
be localized at the bottom of the KS throat, but rather up the throat, at the junction between 
this throat and the ambient Calabi-Yau. Hence, the extra energy brought about by an anti-D3 
brane will contain both a redshifted contribution, coming from the bottom of the KS throat, and 
a non-redshifted contribution, from the junction region. 

It would be interesting to explore whether this extra contribution to the energy could invalidate 
the KKLT construction of de Sitter vacua in string theory [18] . We would like to remind the reader 
that in this construction one adds anti-D3 branes to the bottom of a KS throat in order to lift 
the AdS vacuum to de Sitter. The contribution of these anti-D3 branes to the total energy needs 
to be small in order not to overrun the quantum corrections that are needed to stabilize the 
Kahler moduli. In [18] it was assumed that the contribution of these anti-D3 branes can be made 
exponentially small by placing them in an arbitrarily-long KS throat. Nevertheless, if anti-D3 
branes source non-normalizable modes, the energy they bring cannot be made exponentially small 
by lengthening the throat. Hence, in some circumstances and for certain non-normalizable modes, 
this energy may be too large to allow the Kahler moduli to be stabilized. This without doubt 
negatively impacts the existence of large parts of the landscape of de Sitter vacua in string theory. 

On the other hand, as we will show in Section [3], the force felt by a probe D3 brane in the 
background sourced by anti-D3 branes does not depend at all on the presence of extra non- 
normalizable modes. It is always of the form 1/r 5 , exactly as needed in the KKLMMT model 
for inflation in string theory |13j . 

1.1 Strategy and Results 

Our goal is to identify the modes sourced by anti-D3 branes, smeared in an SU{2) x SU(2) x 
^-invariant way at the tip of the warped deformed conifold. This involves solving the eight 
second-order coupled nonlinear differential equations that govern the metric and fluxes of the 
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Papadopoulous-Tseytlin ansatz [16J. We simplify the problem by solving this system of equations 
perturbatively, treating the anti-D3 branes as a perturbation of the Klebanov-Strassler solution, 
and using the fact that the second-order equations factorize into sixteen first-order linear ones 
out of which form a closed system, that are relatively simpler to solve [17]. One of the main 
achievements of our paper is that we are able to solve these equations exactly in terms of integrals. 

Having set up this rather heavy computational machinery, we then identify which parameters 
correspond to non-normalizable modes and which do not. Our first strategy is to enforce the correct 
UV boundary conditions appropriate for the addition of N anti-D3 branes in the KS background 
and investigate the resulting IR behavior. This amounts to requiring that there are —N units of 
additional Maxwell charge, that there is a non-zero force on a probe D3-brane in the UV and the 
charge is the only non-normalizable mode which we allow for. We find that the resulting infrared 
geometry contains fields that diverge in a subtle manner which seems to be incompatible with 
anti-D3 brane sources. 

As we will discuss in detail in subsection f)6.1.ip these singularities do not appear to have a 
distinct physical origin, as opposed to say singularities that come from explicit D-brane sources. 
On the other hand, they are milder than the latter, as their divergent energy density integrates to 
a finite action. However, even such mild singularities must normally be excluded in order for the 
gauge/gravity duality to make sens^E Whether this singularity is to be allowed or not gives rise to 
two very distinct conclusions. If the singularity is to be allowed, then we argue that we have found 
the correct non-supersymmetric linearized solution for anti-D3 branes in KS. If the singularity is 
to be deemed inadmissible then we argue that there is no solution which can interpolate between 
the KS boundary conditions in the UV and anti-D3 brane boundary conditions in the IR. 

Our second strategy is to enforce anti-D3 brane boundary conditions in the infrared while 
disallowing the aforementioned singularities and examine the consequences for the ultraviolet. We 
show that the ultraviolet of this solution cannot match that of the KS solution with no non- 
normalizable modes turned on. As a consequence we argue that in this case the correct UV 
behavior is a large deformation of the KS boundary conditions. 

There can of course only be one correct solution for anti-D3 branes in KS and since the two 
conclusions we have offered here differ so wildly from each other that it is of utmost importance to 
rigorously settle the issue of admissibility of the singularity which we have uncovered. We present 
various arguments for and against this singularity in the discussion section. 

This paper is organized as follows: In section 2 we set up the perturbation theory following [T7] . 
In section 3 we compute the force on a probe D3-brane in our supergravity ansatz and show that 
is has a simple universal form. In section 4 we solve for the modes and present their perturbative 
expansions. In section 5 we discuss the boundary conditions for D3 and anti-D3 branes in our 
backgrounds. In section 6 we attempt to construct the backreacted solution corresponding to a 
stack of smeared anti-D3 branes in the KS geometry. We conclude and discuss the implications 
of our results in section [7J Appendix |X] presents our conventions and appendix |B] shows the 
asymptotic expansion of the auxiliary fields. 

3 The best example of an unphysical singularity whose action is finite is the negative-mass Schwarzschild solution. 
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2 Perturbations around a supersymmetric solution 



To solve the IIB supergravity equations perturbatively around a supersymmetric solution, we will 
use the method developed by Borokhov and Gubser [T7j . who reduce the set of n second order 
equations for n fields <p a depending on a single (radial) variable r to 2n first order equations for <j) a 
and their "canonical conjugate variables" £ a . For this procedure to work, it is essential that the 
symmetries of the supergravity problem are sufficiently strong so as to just allow a dependance of 
the fields on a single, radial, variable. 

Since we are perturbing around a known solution, the n linear second order equations from 
linearizing Einstein's equation can of course be rewritten as a system of 2n first order equations 
without losing any information []. While this can obviously always be done, the formalism of [17] 
allows eight of these equations to form a closed subsystem when perturbing around a supersym- 
metric solution. Related perturbation theory working directly with the second order equations 
appeared in [19]. 



2.1 The First Order Formalism 

The starting point is the one dimensional Lagrangian 

1 dp dtf 

which we require to have the simplifying property that it can be written in terms of a superpotential 
r l r fdr l rac dW\ ( dr 1 dWx ldW 

where 

1 udW dW 

The fields (f) a are expanded around their supersymmetric background value 0q (which will 
correspond in our case to the Klebanov-Strassler solution [3]) 

P = <f> a + <j>l (X) + 0(X 2 ) , (4) 

where X represents the set of perturbation parameters, and 0° is linear in them. The supersym- 
metric solution (p a = 0q satisfies the gradient flow equations 



dr 2 ^ b 



o^ttd- = , (5) 



while the deviation from the gradient flow equations for the perturbation 0" is measured by the 
conjugate functions £ a , given by 

t s GMa) (M - M W ) , M> d , \°p (G*™) . (6) 

The £ a are linear in the expansion parameters X, hence they are of the same order as the and 
when all the £ a vanish the deformation is supersymmetric. 
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The main point of this construction is that the equations of motion reduce to a set of first order 
linear equations for (£ a ,0 a ): 

0, (7) 
G a % . (8) 

Note that equations (JHJ) are just a rephrasing of the definition of the £ a in (jSJ), while Eqs. © imply 
the equations of motion [IT] . 
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2.2 Papadopoulos-Tseytlin ansatz for the perturbation 

The KS background has SU(2) x SU(2) x Z 2 symmetry. We are interested in a solution for the 
backreaction of a smeared set of anti-D3 branes and we have the liberty to smear these branes 
without breaking the Z 2 symmetry (which exchanges the two copies of SU(2)). Furthermore 
one can see from the supergravity equations of motion that the anti-D3 branes and the fields 
they source do not create a source for the axion Cq. So we are looking for a family of non- 
supersymmetric solutions with SU(2) x SU (2) x Z2 symmetry which are continuously connected 
to the KS solution. The most general ansatz consistent with these symmetries was proposed by 
Papadopoulos and Tseytlin (PT) [16] 

ds \ = e 2A+2 ^ x dsl, + e^dr 2 + e*+\g\ + g 2 2 ) + < ' «(,/ 2 + g 2 4 ) + e" 6 ^ 2 , (9) 

where all functions depend on the variable r. The fluxes and dilaton are 

H 3 = l(k-f)g 5 A(g 1 Ag 3 + g 2 Ag4) + dTA(f'g 1 Ag 2 + k'g 3 Ag 4 ), 
F 3 = F 9l Ag 2 Ag 5 + (2P - F) g 3 A g 4 A g 5 + F'dr A { 9l A g 3 + g 2 A g 4 ) , (10) 
F 5 = J- 5 + *J- 5 , F 5 = (kF + f(2P-F)) gi Ag 2 Ag 3 Ag 4 Ag 5 , 
$ = $( r ), C = 0. 

where P is a constant while /, k and F are functions of r, and a prime denotes a derivative with 
respect to r. There are several different conventions in the literature for the PT ansatz, we compare 
our conventions to those of others in Appendix |A] 

We will denote the set of functions <f) a , a = 1, 8, in the following order 

<p a = {x,y,p,A,f,k,F,$) . (11) 

The metric in the one dimensional Lagrangian ([TJ) is 

G ab <P' a <P' h = e 4p+4A (V 2 + Kj' 2 + 6p' 2 - 6A' 2 + le-*- 2 ^- 2 ^/' 2 + e^k' 2 + 2e 2 *F /2 ) + i$ /2 ^) 

and the superpotential is 

W{4>) = e AA ~ 2p ~ 2x + e iA+Ap cosh y + L 4A+ip ~ 2x (f (2P - F) + kF) . (12) 
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At this point it is worth emphasizing an important point regarding the choice of definition 
for the radial coordinate r. If one redefines r in the ten-dimensional ansatz, or equivalently in 
the Lagrangian (jTJ), this will lead to a significant alteration in the equations of motion (jTJ) and 
(jHJ) since the metric G a b and the connection M a b transform nontrivially. Moreover this also leads 
to an alteration in the definition of the £ a (El). Alternatively one might choose to redefine the 
radial coordinate after having derived the equations of motion (I7j) and (jHJ) in which case the only 
transformation comes from the derivative with respect to the radial coordinate so essentially the 
other terms just get a multiplicative factor of (df/dr) -1 . Of course both such redefinitions are 
completely legitimate, however they are not equivalent since the latter does not alter the definition 
of the £ . 

This is an important consideration since a judicious choice of radial variable can lead to a 
significant simplification in the equations of motion. This must be kept in mind when comparing 
the radial coordinates in |T7] and [20] . The authors of [20J found a particularly natural choice such 
that the equations simplify somewhat and we will employ the same choice in the current work (see 
Appendix lAl for more details). 



2.3 The zero-th order solution: Klebanov-Strassler 

Since we are expanding around the Klebanov-Strassler solution, we tabulate the form of the zeroth- 
order functions here: 

e x ° = i/i 1/2 (|sinh(2r) -r) 1/3 , 

e yo = tanh(r/2) , 

24(±sinh(2r) -r) 1 / 3 



e 6 P0 



e 6A 



h sinh t 



h{\ sinh(2r) - r) 2 / 3 sinh 2 r , (13) 



(r coth t — 1) (cosh r — 1 
Jo = — " 



k = -P 
F = P 



sinh t 

t coth t — 1) (cosh r + 1) 



sinhr 

(sinh r — r) 



sinh r 
% = 

where the function h(r) related to the warp factor is given by the integra@ 

= ho- 32P 2 f tcatht ~ 1 n sinh(2t) _ t y/3 dt n 4 ) 

Jo sinh 2 1 2 

There is an important difference between the Klebanov-Tseytlin (KT) solution [11] and the 
Klebanov-Strassler (KS) one [3J. The KT solution has a warp factor 

27tt (n + aP 2 ln(r /r ) + aP 2 /A 



h KT (r) = ^ — L a = 3/(27r) (15) 



4 We have h = 32P 2 J °° TC ?^I~ 1 {^ sinh(2r) - rf^dr = 18.2373P 2 . 
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which has a naked singularity at some finite distance r s where h becomes zero. The leading UV 
behavior of the KT solution does in fact agree with the KS solution with the identification r = e T ^ 3 
but the subleading behavior of the KT solution is ambiguous due to ability to absorb a shift into 
the singularity For this reason it is not reasonable to compare the subleading behavior of the 
KT and KS solution and in fact deformations of the KT solution (in particular the deformation 
proposed in [12]) which have subleading terms are not in general solutions of ([7]) and (JSj) with the 
0q of the KS solution. 

2.4 £i Equations 

The first step is to solve the system of equations (J7J) for 6- In [17] this system was solved per- 
turbatively in the UV and IR seperately. Since our goal in this work is to connect the UV and 
IR asymptotics, that analysis will not suffice for our purposes. In [20J this system was solved 
analytically with the assumption that 6 = £3 = £4 = 0, which remarkably corresponds exactly to 
the non-supersymmetric deformation generated by a mass term for the gaugino in the dual field 
theory (often referred to as the N = 0* flow). This solution is still not general enough for our 
purposes, we need to lift all of these constraints and consider the most general solution space. 



We find that to solve the ^-equations in general, it is convenient to pass to the basiaj 

6 = (36 - 6 + 6, 6, -36 + 26 - 6, -36 + 6 - 26, 6 + 6, 6 - 6, 6, 6) • (16) 

The equations in the order in which we solve them, are 

g = e - 2 *°(2P/ -F (/o-A;o))6 (17) 

g = -e- 2 *U2P/o- Wo-M)6 (18) 

g = —PeT^ii (19) 

g = -6-^e- 2 *°(P-F )6 (20) 

g = -sinh(2y )6- cosh(2jA))f 6 + Je" 2 * ^ - &o)6 (21) 

g = (Pe 2 ^-sinh(2y )F )6 + (Pe 2 ^-cosh(2y )Fo)6 + ^(/o-A;o)6 (22) 

g = 3e- 2 * - 6 «6 + (5e- 2 *°- 6 "° - e 2 *»(2P/ - P (/o - k )) 6 (23) 
g = 6 cosh sinh y (26 +6 + 6) 

+4 ((Pe 2 *» - cosh(2y )P )6 + (Pe 2w - smh(2y )Po)g) (24) 
where a prime indicates derivative with respect to r. 



5 The inverse is 6 = (±(36 + 6 + | 4 ), 6,6 + 6,-6 - 6, |(6 + 6), ^fe + 6), 6, Is) 
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2.5 (j) 1 Equations 

Once the equations are solved, one can insert the solutions in the $ equations, which we write 



here explicitly. We start by performing a linear field redefinition^) 

a = (x - 2p - 5A, y, x + 3p, x -2p- 2 A, f, k, F, $) . (25) 

where all the fields on the right hand side are the first order perturbations. 
Then the system of equations is (in the order in which we will solve them) 

0' 8 = -4e- 4{Ao+po) is (26) 
0' 2 = -cosh|/o02-2e- 4(Ao+po) e2 (27) 

0' 3 = _3e-^- 2 ^ 3 -sinhyo h - V 4 ^^ (9& + 5f 3 + 2f 4 ) (28) 

o 

$ = 2e- 6 ^- 2 -V3-sinhy o 2 + ie- 4 ^ +p °)(fi + 3f4) (29) 

o 

0' 5 = e 2w (F o -2P)(20 2 + 8 ) + e 2w 7 -2 e -^^ (30) 

= e- 2 ^(F o (20 2 -0 8 )-0 7 )-2e- 4 ^ + ^ +2 ^-^H (31) 

= \(k-k + {k -fS,)-2e-^ + ^ +2x °i 7 (32) 

04 = V 2 ^(/ o (2P-F o ) + A;oFo)(20 1 -203-50 4 ) + ie- 2:Eo (2P-F o )0 5 

+ ^- 2 ^F o 06 + ^-^(fco - h)h - ^e- 4 ^+^)6 (33) 



2.6 Imaginary (anti)-Self Duality of The Three Form Flux 

Before moving on to solving the first order system of equations, it is worthwhile considering the 
conditions for self duality of the three form flux. The complex three-form flux of the KS solution 
is well known to be imaginary self dual (as are all warped Calabi-Yau solutions of IIB supergravity 
[2Tj 122] ) and one might be tempted to consider the possibility that the flux turned on by an anti-D3 
brane is imaginary anti-self dual. It is important to note that since the metric is deformed at first 
order (as opposed to deformations of AdS§ [23]), it is not consistent to consider the self duality of 
the flux with respect to the background metric. Instead, one should analyze the condition for self 
duality using the perturbed metric. 

In our conventions, the condition for imaginary self duality is 

ISD : e $ * F 3 = -H 3 , (34) 

and the condition for imaginary anti self duality is 

anti - ISD : * F 3 = +H 3 . (35) 

6 The inverse is <j) a - </>g = (±(-2</>i + 2<f> 3 + 50 4 ), 4>2, je(^4>i + 3<fe - 504), \{~<t>i + <M, 4>5,4>6,4>7, <t>$) and to avoid 
confusion we abuse notation at from now on put all indices down. 
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In terms of the PT ansatz, these conditions lead to 



e*(2P-F)e 2y = Tf' (36) 

F' = *\e*{k-f). 

with the upper sign for the ISD condition. 

Using Eqs (I30ti26[) and the fact that the zeroth order three-form flux is imaginary self dual, the 
conditions that the first order flux be ISD become 

fe-fe) = 0, 

+ = 0, (37) 
h = o. 

It is easy to see that these immediately lead to £i = and force (£4, £g) to be constant. Hence, the 
generic solution to the system of first-order equations [XT] does not have purely ISD flux, contrary 
to recent claims in [21]. One can also see this by examining the small-r expansion of the dilaton. 
ISD fluxes do not source the dilaton. Allowing for a delta-function source at the origin, the dilaton 
in an ISD background should be equal to 

r 3 



c >+< + S-^+-) <* 



where the constant C2 multiples the free Green's function of the deformed conifold [21]. It turns out 
that the full expansion of the dilaton near the origin also contains even powers of r. Setting those 
to zero and demanding the coefficient of the odd terms to be proportional to those in the Green's 
function (135]) implies certain conditions on the integration constants (Xi = 0, X 5 = — X 6 = X?), 
which are actually weaker than (13T|) . 

It is interesting to note that all BPS solutions of (jTfS]) have ISD three-form fluxes since for 
supersymmetric solutions all the £'s vanish. If one relaxes the Z 2 symmetry in the PT ansatz then 
there exist BPS solutions with non-ISD flux [25] [26]. 

The conditions for the first order deformation to be anti-ISD are of course more involved: 

(£5-6) = e 4 ( Ao+ ^- 2 ^(F o (20 2 -0 8 ) + 7 ), 

(S + S) = e 4 ( Ao+ ^- 2 ^((F o -2P)(20 2 + 8 ) + 7 ), (39) 
ir = -l^ Ao+Po) - 2x °(^-^)-^( k o-fo))- 



One can try arguing that these conditions are necessary in order for the first-order perturbed 
solution to correspond to anti-D3 branes. However, as we will find in the next section there exists 
a much simpler and more physical signature of the presence of anti-D3 branes: the force on a probe 
D3 brane. 
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3 The Force on a Probe D3 Brane 



Before solving the above equations, we compute the force on a probe D3 brane in the perturbed 
solutions. As is well known, such a brane feels no force in the unperturbed Klebanov-Strassler 
solution. It appears at first glance that both the Dirac-Born-Infeld (DBI) and the Wess-Zumino 
(WZ) actions are rather complicated in the perturbed solution, however we show below that by 
using the first order equations of motion, most terms in the total action cancel and the final 
expression for the force on a probe D3 brane is quite simple. 
The DBI action for the probe D3 brane is 

V DBI = V-9009lig22g33 

= e 4A+4p - 2x , (40) 

and in the first-order expansion, the derivative of this action with respect to r is 

dV DBI 



F DBI 



dr 

= _e 4Ao + 4 f«- 2 ^ (1 + AA l + 4pi - 2xi ) (4A' + 4p' - 2x' + AA[ + Ap[ - 2x[) 

= - e 4A ^- 2 ^(-2M^' + 4p' - 2*' ) - 20' 4 ) , (41) 

where V^ 31 is the DBI action in the unperturbed solution, 

F° BI = _ e 4 ^+4po-2xo ( 4y ^ + Ap ' o _ 2x ' Q ) } (42) 

and we have used the definition of 4 in l$Fo§. 

To compute the form of the WZ action one needs to know the value of the RR potential 
along the worldvolume of the brane, which we obtain by integrating the RR field strength F^ 5 \ It 
is however easier to compute instead the derivative of this action with respect to r, which gives 
the force exerted on the D3 brane by the RR field: 

dv wz 
pwz av 

dr 

— - r 0123r 

= —{kF + f(2P - j P)) e 4A+4 P -4x_ 

The zero-th order and first order WZ forces are 

F^ z = -(fc F + /o(2P-F ))e 4Ao+4p °- 4 ^ (43) 

F WZ = _ e 4A 0+ 4 P o-4,op oFo + /o(2 p_ Fo))( 4 Ai+4pi _4 Xi)] 

_ e 4A 0+ 4 PO -4xo [hF() + Fi(A . o _ /o) + /i(2 p _ Fq) ] (44) 

It is not hard to check that the zeroth-order WZ and DBI contributions to the force cancel, as 
expected for a BPS D3 brane in the KS background. 
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Using this, as well as the equation of motion for <j)' 4 we get 

F° BI = _ e 4Ao+4 Po -2*o (-20 4 (4A' O + V - 24) - 2& 
= 2e 4 " 4o+4po_4x ° 



~{k F + /o(2P - F o ))(20! - 203 - 1O0 4 ) 
5 



+e 4A 0+ 4 P0 -4*o [fciFo + Fi(A . Q _ /o) + h( 2p _ Fq)] + ^e- 2 ^ . (45) 



From the definition of the (ft 1 , we have that 

l 4 = 10{A 1 +p 1 -x 1 ), 



thus the first and the second term of the Wess-Zumino force (1441) cancel against corresponding 
terms from the DBI force (j4"5j) . Hence to first order in perturbation theory, the only contribution 
to the force felt by a probe D3 brane comes from the term proportional to £i: 

F = F DBI + F wz = -e- 2 ^^. (46) 

3 

As we will show in detail in section HJ it is rather straightforward to compute the UV asymptotic 
expansion of £i by expanding the integrand in (153"]) for large values of r: 

= Xi3(l - 4r)e- 4r/3 + C( e - 10r/3 ) (47) 

By expanding in the ultraviolet e~ 2x ° as well 

e~ 2xo = — — + 0(e- 2r ) , (48) 

3P 2 (4r - 1) v ; ' v ' 

we can see that the UV expansion of the force felt by a probe D3 brane in the first-order perturbed 
solution is always 

F T ~ X ie - 4r/3 + C(e- 10T/3 ) . (49) 



Recalling that in the UV, r is related to the canonical radial coordinate 



,r/3 



(50) 



it follows that 



and thus the potential goes like 



*~fr + o(;lr] ( 51 ) 

V ~ C + ^ . (52) 

Having obtained this simple universal result, we pause to consider its physical implications. First, 
the force is largely independent of which modes perturb the KS solution and on whether these 
modes are normalizable or non-normalizable; it may be zero in certain solutions (such as [20] . 
where X\ = 0) but when non-zero it has a universal form (T5TT) . Out of the 14 physically- relevant 
SU(2) x SU (2) x ^-preserving perturbation parameters, only one enters in the force. 
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Second, this result agrees with that obtained in [13J. In that paper, the force felt by a probe 
anti-D3 brane in KS with D3 branes at the bottom was shown to scale like 1/r 5 and to be linear 
in the D3 charge. It was also argued using Newton's third law that the same force should be felt 
by a probe D3-brane in KS with anti-D3 branes at the bottom. Our result implies therefore that 
in order to describe anti-D3 branes in the infrared, one must have a nontrivial £1, and furthermore 
that the constant X\ must be proportional to the number of anti-D3 branes. 

We would also like to note that in the ultraviolet the mode (in the fa) proportional to X\ is 
a normalizable mode decaying as r -8 , in fact this is the most convergent of the 14 modes. The 
necessity for having such a mode in order to find the r~ 5 force predicted by [13] was discussed 
in [12] , although the mode itself was not identified. However, the force analysis alone in no way 
supports or disfavors the possibility that an anti-D3 brane sources a non-normalizable mode. The 
force has exactly the same r-dependence regardless of whether non-normalizable modes are turned 
on or not. Hence, the cancelation of terms in the force up to order 1/r 4 is universal. 

4 The Space of Solutions 

In this section we find the generic solution to the system (I17H33I) . which depends on sixteen inte- 
gration constants (although only fourteen are physical). As we will see, the full solution involves 
integrals that cannot be done explicitly but these integrals can be expanded both in the UV and 
IR, giving rise to sixteen UV and sixteen IR integration constants. Having the expression for the 
full solution permits one to relate numerically the UV and IR integration constants. 
The first equation ( 1TTI) is solved by 

£x = X x exp ( jT dr'e- 2xo [2Pf - F (f - k )}^ (53) 

and this integral cannot be done explicitly. In a cruel twist, having argued above that the inte- 
gration constant X\ is the only mode which contributes to the force on a D3 brane and is thus 
the signature of the anti-D3 brane background, we have immediately found that as long as X\ is 
nonzero there is no analytic expression for any of the £i modes. On the other hand, when X\ = 0, 
we can fully solve the system of equations for the £; analytically. We proceed by initially setting 
X\ to zero, finding an analytic expression for all the other £j, and then building on top of this, the 
integral expressions for the full solution with nonzero X\. 

4.1 Solving the ^ Equations for X\ = 

Here we solve the £j equations for the choice 

X 1 = |i = (54) 

and later we will find integral expressions for the full solution including £i. Having set £i = 0, we 
then immediately read off that (£4, £5) are constant 

£4 = X 4 , £ 5 = X 5 . (55) 
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We now have a pair of coupled o.d.e's for (£6, £7) 

§ = ~ir, (56) 
g = -sinh(2 2 / )X 5 -cosh(2 Z / )e 6 . (57) 

Doing the derivative of (156"|) and using (j57p one easily determines the solution 

(1 + 2r)X 5 + X 6 - 2tX 7 + X 7 sinh(2r) 

^ 6 = 2iinh7 ' (58) 

~ ((1 + 2r)X 5 + X 6 - 2rX 7 )coshr - (2X 5 + X 7 (cosh 2r - 3)) sinh r 

2 sinh 2 r 

To ease notation the two integration constants appearing in the second order equation for have 
been denoted X 6 and X 7 , but one should keep in mind that they both appear in £ 6 and £7. To 
compare to the solution in [20J one sets X 7 = — X 6 = X 5 . 
We then directly integrate (l2"2"j) and (1251) to find 

is = (l-e 2r ) _3 (2e 3r (P((-l + 2r + 4r 2 )X 5 + X 7 + 2r(X 6 -2rX 7 ))coshr 

+P(X 5 - X 7 ) cosh 3r + (3X 8 + P(X 5 (1 - 5r) - 2X 6 + (-3 + 5r)X 7 )) sinhr 
-(X 8 + P{1 + r)(X 5 - X 7 )) sinh 3r)Y (60) 



£ 3 = 2X3 ( sinh 2r - 2r) . (61) 
Finally we insert the result for £3 into ([241 and solve the o.d.e. to find 

£ 2 = - ((2X4 + 4r(-4X 3 + 3PX 7 )) cosh r + 2(6X 2 - X 4 ) sinh r) 

+77^7^ ((2(^5 - X 7 ) - (X 5 (l + 4r) + X 6 - 4rX 7 ) 
2 sinh r 

+ o PT „ (X 5 (l + 2r) + X 6 - 2rX 7 ). (62) 
2 sinh r 

So far we have introduced 8 integration constants 

(X 1 ,X2,X 3 ,X4,X 5 ,X 6 ,X 7 ,X 8 ) (63) 

however the the zero-energy condition is a constraint amongst these 

- X 4 + 6X 2 - 3PX 5 - 4X 3 + 9PX 7 = . (64) 

This leaves 7 independent constants, but for the time being we will keep the 8 constants since this 
will provide a clearer picture of the UV asymptotics. 

4.2 Solution for X\ 7^ 0, UV and IR expansions 

When £i is nonzero, one can find integral expressions for the solution, but the integrals cannot 
be performed explicitly since £1 itself, given in (153|) . cannot be found explicitly. £1 appears on 
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the right hand side in all of the £ equations (except for the equation for £§, but the right hand 
side of ( 1221) contains £5,6,7 which are modified by £1) and therefore all £'s get an extra contribution 
proportional to X\. A simple way to understand this extra contribution, is to recall the following 
method to solve a system of coupled o.d.e's, which will in fact be used to solve the equations for 
0. Given the system 

g^A/g. + h (65) 

where i = l,...,n, and g,A,b are functions of a single variable r, one first finds the solutions 
g 3 hi of the homogeneous equations (again j = 1, . . . ,n). Then one constructs the ansatz for the 
inhomogeneous solution by a linear combination of the homogeneous solutions where the coefficients 
are promoted to functions, and then it is easy to see that one has 



(66) 



The integration constant here corresponds to adding any amount of the homogeneous solution to 
g>j. For the simple case n = 1, that will appear repeatedly in the following, we have 

g(r) = X(r)g h (r) , g h = e! A , A = f ^\dr' , (67) 

Jr 9h{r) 

where the integration constant is related to the choice of tq. 

Adding £1 on the right hand side of (118ti24p amounts to the change bi — > bi+Xifi(r) = b^ + Xibj 
where b\ contains £1 0. This implies 

\ \ + X\ [ T dr'ig^b) = \ Q i +X 1 X\ , g i ^g i + X 1 g J h i X] . (68) 
Jo 

Note that in the expression for A we have chosen to absorb the n overall integration constants in 
A° and fixed the limits of the integrand in Xj. 

The integrands can always be expanded around the IR and UV, and the integrals can be 
performed in these limits. Expanding ( 153]) for small r amounts to expanding the integrand in a 
series of powers of r, doing the integral explicitly and then expanding the exponential agairH This 
gives 



~tr tr( 16(f) 1/3 P 2 t 2 

Zi R = x[ f 1 ^_ + (, 



(69) 



The contribution from £1 on the rest of the £'s can be obtained similarly by expanding the inte- 
grands in A". We give the IR behavior of all the £'s in Appendix lB.il 

The UV behavior is slightly more tricky. On one hand, the expansion is in powers of e~ T ^ 3 ~ 1/r, 
where at each power there are terms proportional to r n ~ (logr) n (typically only n = and n — 1 



7 The indices and 1 in b are used to distinguish between the expression without £1 and that purely due to £1 . 
They should not be confused with zero and first order in perturbation theory (they are both first order). 

8 In general, if for a given expression to be integrated in the IR there are negative powers of r showing up 
in the expansion, we subtract the infinite contribution from the lower limit of the integral, namely we replace 

Jo" fnt- n -> f fnt- n = (for U £ 1). 
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appear). Fortunately all these terms can easily be integrated. On the other hand, expressions like 
( )53l) are split into 



POO 

6 = X ie xp I J dr'f(r') j = exp (/ dr' fir') j exp U dr'fir' 

ee X^ v ex P ( [ T dr'fir')) (70) 



where the integrand can now be expanded around large r, and the integral can be done analytically. 
The integral from zero to infinity giving the relation between X± and X\ = X[ can in principle 
be done numerically. This gives 

= X tAV e -4r/3 (3 _ 12r) + O ( e -10r/ 3) (n) 

For the other £'s, one uses the same expansions to obtain the UV behavior of Xj. We get from 

AMD 

X x X\ = X x QT dr'igj^b) + | T dr'fo; V &}) = *i + / T ^'(^ V&i) , (72) 

where fcj are some constants that can be obtained numerically, and in the remaining integral one 
expands the integrand for large r and performs the integral analytically. The contributions X\ki 
shift the constants of integration in A°. We therefore define 

X™ = Xl* + Xika . (73) 
We give the UV expansion of all the £'s in Appendix IB. 21 

4.3 Solving the (jf Equations 
4.3.1 The Space of Solutions 

Here we show how to solve the system of <p l equations (I26H33I) . assuming X\ = 0. First we note 
that <f)$ is given by the following integral, which however cannot be solved analytically 

^ = - M / dT (|sinh(2t 8 )-^/3 +y " (74) 

where £g is given in ( 160]) . 

The equation for <fi 2 is solved using (1671) . which implies 

A 2 (t) ~ 9 f , sinhrfo rf? . , 

smhr 7 (| smh(2r) — r) 2 ' d 

where (sinhr) -1 is the solution to the homogeneous equation, and £2 is given in (1621) . The same 
trick is used for 3 , which is given by 

(76) 



sinh(2r) - 2r 



. , .02 16 sinh 2r - 2r 5X 3 + X 4 \ , B 

A 3 = / dr sinh 2r - 2r — ^ ^ — ^- / f 0/Q + K, /iJ 

' v ; M sinhr 3 ± sinh 2r - r 2 /3 ' 
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We obtain an integral expression for 



a. 



dr 



8 sinh tc 



+ 



+ 



32 



3 sinh(2r) - 2r sinh r 9 (± sinh(2r) - r) 2 / 3 
and then using ( |66|) . we find that the fluxes (05, 06, 07) = (/, k, F) are given by 



IR 



(77) 




/2r - sinh r - tanh(r/2) - 
2r + sinhr - coth(r/2) + 



i 



2cosh 2 (r/2) 2cosh 2 (r/2) 

T 1 



V 



sinh -, 



— cosh T 



2sinh 2 (-r/2) 



2sinh 2 (r/2) 
1 



sinh i 



1 




(78) 



This 3x3 matrix contains the solutions to the homogeneous equations, and the derivatives of the 
functions A a are given by 




4 sinh t 

• | (cosh T — 



sinh i 



cosh T 

2 sinh 2 r 
(sinh(3r)+4r cosh r— 7 sinh r) 



sinh i 



sinh - 



8 sinh 2 i 

:-i+r 



cosh T N 
sinhr ' 




(79) 



where 65, b^, 67 are the rhs of (|30|) . (|3T|) and (|32l) respectively, and the 3x3 matrix is the inverse of 
the one in (I7HI) . We will call the constants of integration in these equations Y 5 , Y e , Y 7 , even though 
the three functions depend on the three of them (i.e. the integration constant in 05 is not just I5, 
but a combination of 15,6,7)- 

Finally, the equations for 4 is solved using the same trick, namely (|6TI) except that the solution 
to the homogeneous equation can be found only as an integral expression. Formally, it is 



04 = A40, 



1/f. 



A, 



— + Y IR 

Hh 



30) 



where 64 is the right hand side of (1331 (setting 4 to zero) and 

rsinh(3r) — cosh(3r) + 5r sinhr — coshr(4r 2 — 1) 



Ml, 



cxp 



8P Z / dr 



h(r) sinh r(cosh r sinh r — r) 2 / 3 



(81) 



4.3.2 IR behavior 



We now give the IR expansions of the 0\ We write the divergent terms and the constant terms 
since terms which are regular in the IR will not provide any constraints on our solution space. 
Here we have used the zero energy condition (|84|) to substitute for 
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1 

^ (2/3) 1/3 (6X 8 /fl + P{5X' R - Xl R - 6X' R )) + Y 8 IR + 0{r) 
, lo S r 



-V/*4 
r t 



l6(2/3) 1/3 (-(l + l6(2/3) 1/3 ^)Xl R +(-2Xl R +PXi R +PXi R )) + 0(r) 

V 3 iIq J 



t 13Y 3 IR l f Y 2 IR 3Y 3 IR X( R QA(2/3) 2 / 3 P 2 



r3 4 + 7v^ 20" + h P{X " +X& H(2/3) v 

f^ll (-Xf* (8(2/3) 1/3 /3 + (2/3) 2/3 128P 2 //i ) + 8(2/3) 1/3 (-2X 4 /K + P(X IR + X 6 ))) + 0(r) 
1Y* R 1/ „, rR V, /fi ,„/ „„ 2 2 / 3 3 1 / 3 128P 2 



3 ^_ + i( - 2Y/* + ^- + X^(l0(2/3) 4 / 3 , 



+4(2/3) 1 / 3 (10Xf - 6P(Xf + X™))) + (x i /R (l6(2/3) 4 / 3 + ( 2 / 3 )^ 512p2 
+X 4 /K 64(2/3) 1/3 - (X 5 /R + X 6 /i? 32(2/3) 1/3 )) + + O(r) 



^ + Y 7 IR + X IR ^f ( 16(2/ ^ )2/3jP2 - (2/3) 1 / 3 ) - Xf 16(2/3) 1 / 3 P + C(t) 



'o 



J_ ( _ 2 r 6 /i? + (2/3) 1 / 3 32PX 1 /i? - (4/3) + X 6 /i? 

+ (*Jl! + y/« + Xf (2/3) 1 / 3 16P - + y (1/3) 1 / 3 P 2 )X 5 /R - 32(2/3) 1 / 3 PX 7fl 

iy IR - f— + — 
3 2 V 9 3 



1P y ir + ^(i/ 3 ) 1 /3p2j x /i ? + 32 (2/3) 1 / 3 P 2 (X 6 /i? /3 + X. 



IR-* 
7 t 



(r R A2/3) l ' 3 32P 512(2/3) 2 / 3 P 3 32(2/3) 1 / 3 P fn WJ7 „ 2/WJ? w ^ 
+ logrfX^^^ + LLJ ) + (2X IR - P 2 {Xi R + X' R )) 

+0(r) 

1 f ^ir , h 



T 



( - Y'« + |(A' 5 ' R + + 0(t) 



I , 16(2/3^/* 2 V3 3 V3 8Pr/fi + Y , s(8(2/3)I/3 + 64(2/3)2/3p2/fto _ 4096^ 



+X , R 256(3/2) 1 /3 P 2 + 8(2/3)4/3(x ^ + x , R) _ 16(2/3)V3p2 ^ + ^ + ^ 



4.3.3 UV behavior 

Here we provide the UV asymptotics for the (fi l including terms which fall off as e _4T//3 and slower. 
However as shown in the table I4.3.3[ certain modes have leading behavior in the UV which is more 
convergent than this. 
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8 = Y 8 UV + 2 1 / 3 24e- 4r / 3 (P(3 + 4r)(Xf y - X^) + 4X 8 c/y ) + 0(e' 10 ^ 3 ) 

2 = -2 1 / 3 48e-^ 3 (2X 2 7y - (-3 + 2r) {2X^/3 - PX^ + 2e~% uv + 0(e' 7T / 3 ) 

3 =-2 1 / 3 10X 3 7y e 2 ^ 3 - e^ 3 2 10 / 3 (36X^ v + 2Xf y + X 3 VV '(112 -^r) - 36PXf y (3 - r)) + 0(e"^) 

^ = _ 2V320X3 7 V^ 3 + 2^4 (lOSX^ + Xf y + 176X3^ - lSQPX^) e" 4 ^ 3 

_ e -4r/3 r2 i/3 16 ^ x t/y _ 27 PX 7 /y ) + 0(e- 2r ) 
vuv 

05 = -\e* + {-Y™ + Y™) + - PY™)t 

+2 1 ' 3 ^{2PX^ V + 7PX 3 l/y + 7P 2 Xr)e- r/3 - lil^L^V 
+^5>f y + AY 6 UV + 2PY 2 UV - 2PY 8 uv y~ T + i( - AY 5 UV + 2PY 8 f7y ) e -V 
+2 1 / 3 72P( - AX™ + X 8 C7y - 5X™ + 2P(X 5 l/y + X 6 c/y ))e- 4 ^ 3 

-2 1 / 3 72P(4X 2 7y - 2PX" v - + PX!/ v y~^ 3 r 

-2 1 / 3 72P( - + 4PX 7 /y )e- 4 ^ 3 r 2 + 0{e- 2T ) 

vuv 

06 = +^-e T + (-Y 5 uv + Y 7 uv ) + (2Y 5 uv -PY 8 uv )r 

-2^36 {2PX™ + 7PX 3 vy + 7P 2 X 7 /y )e-^ 3 + Ili^L^/ 3 r 

_I + 4Y uv + ^t/v _ 2 pr^v j e -r _ I ( _ 4y uv + 2 py 8 ^ e -r r 

+2V 3 72P( _ + X^ - 5X™ + 2P(X 5 l/y + X^e'^ 3 

-2 1 / 3 72P(4X 2 c/y - 2PXf y - + Pl^je-^r 

-2 X / 3 72P( - + 4PX 7 /y )e- 4T / 3 r 2 + 0(e- 2T ) 

Y uv / \ 

7 = |-e r + 2 1 / 3 36P^-6X 2 7y -(9-4r)X 3 7V ' + (10-4r)PX 7 7y Je- r/3 

+ (i(4r - 1)Y 5 UV - 2Y 6 UV + PY 2 + P(r - l)>f y )e- r + 0(e^ /3 ) 
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04 - 3(4r-l) 6 + (4r-l) 6 

-AY 7 uv + (3 + At)Y 8 uv l_ uv _ 2Y[^ 32Y4 UV (12 - 85r + 25r 2 ) _ 2r/3 
2P(4r-l) P 5 5 1125(1 - 4r) 2 6 

nl/3 / 

+ ¥ — IV ~ 18P ( n + 8r ) X 5 /l/ + 72 ( X + 8t ) X 2 V + 2 (-5 + 8t)X" v - 72X% V 

+9P(95 - 56, + 80rW - 6«3<-3279 + 17012, - 20785r* + 14900r')y 4r/3 + _ 2 

375(4r — 1) / 

To understand the holographic physics of the ej) 1 modes it is useful to tabulate the leading UV 
behavior coming from each mode. For each local operator Oi of quantum dimension A in the 
field theory, the well known holographic dictionary [271 12H] relates two modes in dual AdS space, 
one normalizable and one non-normalizable. These two gravity modes are dual respectively to the 
vacuum expectation value (VEV) (0|Cj|0) and the deformation of the action 5S ~ f d d xOi. 

normalizable modes ~ r~ A -h- field theory VEV's 
non — normalizable modes ~ r A ~ 4 <H- field theory deformations of the action, 

and we recall that we have in the UV r = e T ^ 3 . For backgrounds like Klebanov-Strassler which 
are asymptotically AdS only up to certain logarithm terms, it is expected that this dictionary still 
holds. In Table 1 we have summarized which integration constants correspond to normalizable 
and non-normalizable modes. It is very interesting to note that in all cases a normalizable/non- 
normalizable pair consists of one BPS mode and one non-BPS mode. 



Table 1: The UV behavior of all sixteen modes 

in the SU(2) x SU(2) x Z 2 symmetric 
deformation of the Klebanov-Strassler solution. 



dim A 


non- norm / norm 


int. constant 


8 


y»4 / ry 8 


Y A /X l 


7 


ty^ j ry f 


Y 5 /X e 


6 


ry^ 1 ry '6 


X 3 /Y 3 


5 


ry j ' ry 5 




4 


yiO / ry 4 


^8) ^1/ -^5, -^4, V 8 


3 


ry 1 / ry 3 


X2, X 7 /Y 6 , y 2 


2 


y 2 /y 2 





It perhaps bears repeating here that the are integration constants for the £j modes and 
break supersymmetry, while the Yj are integration constants for the modes </>\ One key result from 
this table which cannot be gleaned from the field expansions we have provided, is that the mode 
£1, whose integration constant is Xi and which is the only mode responsible for the force on a 
probe D3-brane, is the most convergent mode in the UV. 
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4.4 The Zero Energy Condition and Gauge Freedom 



In addition to the first-order equations of motion (jSJ) there is a constraint, coming from the R TT 
component of Einstein's equation: 

i drd<p b _ 

2° ab ^^ ~ V{<P) (82) 

which must be enforced in order for the solutions of our system of equations to correspond to 
solutions of supergravity. This constraint is imposed after fixing the gauge freedom associated 
with redefining the radial coordinate. When linearizing the PT ansatz around the KS background 
this becomes a constraint on just the £ a 

6 £-a (S3) 

In addition, the common rescaling of the coordinates (t, x%, X2, £3) can be used to eliminate the 
integration constant in the metric function A (since A = |(— 0i + 04 ) in the UV this integration 
constant is Y^ v ). Hence the linearized solution depends on an even number of physically- relevant 
parameters (fourteen), consistent with the expectation from the AdS-CFT correspondence. 

On our solution to ([7]) and (jSJ), we find that the zero energy condition is a linear algebraic 
relation amongst the integration constants 

- 6X 2 + 4X 3 + X 4 + 3X 5 - 9X 7 = (84) 

(note that X\ does not enter in the condition). Using this to eliminate X 4 or X 5 then leaves just 
two normalizable and two non-normalizable modes with A = 4 in Table 1. 



5 Boundary Conditions for Three-Branes 

Within the space of solutions we have derived in Section H] we are interested in finding the modes 
which might result from the backreaction of a collection of smeared anti-D3 branes. These branes 
are placed at r = and smeared on the finite size S 3 , and for describing them it is necessary to 
carefully impose the correct infrared boundary conditions. 

The gravity solution for a stack of localized D3-branes in flat space has a warp factor h(r) = 
1 + Q/r 4 and as r — > the full solution is smooth due to the infinite throat. However when these 
branes are smeared in n-dimensions, the warp factor will scale as r~ 4+n as r — > since it is now 
the solution to a wave equation in dimension d = 6 — n. This is the IR boundary condition we will 
impose on the metric. 

In addition we must impose boundary conditions on the various fluxes. This is relatively 
straightforward for D3-branes in flat space, where the energy from cancels against that from 
the curvature. In the presence of other types of flux, the IR boundary conditions are slightly 
more subtle. When the background is on-shell, contributions to the stress tensor from all types 
of flux taken together cancel the energy from the curvature: this is the basic nature of Einstein's 
equation but this is too loose a criterion to signal the presence of D3 branes. The correct boundary 
conditions for D3-branes should require the dominant contribution to the stress-energy tensor to 
come from the F^ flux. 
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5.1 BPS Three-Branes 



The Klebanov-Strassler background is a smooth solution and has no D3 Page charge. The D3 and 
D5 Page charges are 

(fig* = L_ f (jr(5)_ jB (2) Aj p(3) + l S (2) AjB (2) AF ( 1 ) 

(4tt ) Jm 5 ^ 2 

Q%*° = ^ F (3)_ jB (2) AF (l) 



I S3 

and evaluated on the KS background these are 

Qd3 = 0, Q D5 = P. 

It is trivial to obtain the BPS solution in which smeared D3 branes are added at the tip of the 
deformed conifold. In the PT ansatz one can add D3-brane sources by shifting J 7 ^ by a constant 
multiple of the volume form on T 1 ' 1 . This shifts the D3 Page charge, as well as *J r( - 5 \ in such a 
way that the warp factor becomes singular at the tijH 

The integral that gives the warp factor remains of the same form as in the original KS solution, 

( r f„(2P - jo) + k F , 
" (T W smhVA'V) <85) 

with K{r) = ^ 2i/ 2 3 T s 7nh"r ' ' Shifting J 7 ^ by a constant corresponds to shifting f and k by equal 
amounts. Under a shift of N/ (2P) the D3 Page charge shifts by 

Qds -> Qds + N (86) 

while Qob remains invariant. This introduces in the IR a 1/r singularity in h(r), whose physical 
interpretation is obvious: we have smeared BPS D3-branes (whose harmonic function diverges as 
1/r 4 near the sources) on the S 3 of the warped deformed conifold. 

In order to "calibrate" our perturbation-theory machinery it is interesting to see how this BPS 
solution can be obtained in an expansion around the BPS KS background. First, all the fields, 
that correspond to modes that break supersymmetry at first order, should be set to zero, and hence 
Xi = 0. Furthermore, since £j are zero, Y^ R = Y^ v , and requiring no divergent terms in the IR of 
02 sets these to zero. The same thing happens with Y 3 and Y\ (for the latter we require no constant 
term in the UV of 0i). The constant Y 8 has a very clear meaning: it corresponds to a shift in the 
dilaton. Note that this shifts 5 , 6 and 4 as expected from the exact KS solution^. We consider 
perturbation where the dilaton is not shifted at infinity and hence set Y% R = Y% = 0. The fields 
05,06 and 07 then satisfy the homogeneous equation, and their solution is obtained from flTHj) by 
setting A5 ; 6,7 to ^5,6,7- Requiring no exponentially divergent terms in the UV and no divergencies 
in the IR set respectively Y 5 = and Y 6 = 0. The only non-zero integration constant one can have 
is therefore Y7. 



9 Note that one can also shift the D3 Page charge by an integer multiple of P, by a large gauge transformation 
of i?2- However, this will not affect *J 7 ^ or the warp factor, and hence will leave the physics invariant. 

10 The NSNS B-field has a factor of g s which we have set for simplicity to 1 in Eq. (|13|) . Similarly the integral 
that gives h in Eq. (fT4|) contains a factor of g^; shifting g s in the BPS solution changes <j>4 by a term proportional 
to t/(4t — 1), exactly as in our UV expansions. 
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Hence, the perturbation corresponding to adding N BPS D3 branes is obtained by just setting 
Y? ~ N. Note that this perturbation causes the warp factor 04 to diverge in the infrared as N/r, 
exactly as one expects from the full BPS solution. All the other fa change by subleading terms, 
except 05 and 06, which as argued above shift by N. 

Note that in the UV Yj multiplies a non-normalizable mode corresponding to the fact that we 
have changed the rank of the gauge group of the dual field theory by adding D3-branes: 

= -3P(4r-l) + - 

Hence the new warp factor is 

h = e 4A+4p-2* = pV 4r/3 (4r-l)(l-20 4 ) 

= PV 4 ^ 3 (4r - 1) + l^e" 4 ^ 3 + ... , 
and we can see that Y 7 multiples a 1/r 4 term, as one expects for pure (non-fractional) D3 branes. 

5.2 Anti-BPS 3-branes in Anti-Klebanov-Strassler 

The anti-KS solution is obtained from the KS solution by flipping the sign of and J 7 ^ but 
not of F^ 3 '. This flips the signs of ko and fo in the expression for and has the same D5 
brane Page charge as the KS solution, namely Q^ 9e = P- It is rather straightforward to see 
that this solution is also supersymmetric but it preserves different Killing spinors than the original 
KS solution. With the addition of iV anti-D3 branes smeared on the S 3 the warp factor h again 
diverges in the infrared as — , and the functions / and k are equal and approach — N. The D3 

brane page charge is Q^ e = —N. 

The KS solution with D3 branes at the bottom can be thought of as coming from first taking 
the Ricci-flat deformed conifold with a topological flux, puting D3 branes in, and backreacting 
this configuration. The orientation of J 7 ^ (and subsequently that of H^) is determined by the 
orientation of these D3 branes. In a similar fashion, one can obtain the anti-KS solution with 
anti-D3 branes by starting from the same Ricci-flat deformed conifold with the same topological 
flux and putting anti-D3 branes in; this backreacted system has J-( 5 ' and of opposite 
orientation compared to the solution in section 15.11 Hence, for one orientation one obtains BPS 
D3 branes in the KS solution, while for the opposite orientation one obtains the anti-BPS 3-branes 
in anti-KS. 

Note that if one thinks about the solutions in this way, it does not make sense to a-priori fix 
the sign of the charge of the "fractional" D3 branes. The only quantity whose orientation is fixed 
is the integral of F^ on the 3-cycle. Hence, from this perspective, changing the orientation of 
the D3 brane sources reverses the sign of JF^ throughout the solution, and hence transforms the 
"fractional" D3 branes into "fractional" anti-D3 branes. 

5.3 Anti-BPS 3-branes in Klebanov-Strassler 

If one now inserts anti-D3 branes in the BPS Klebanov-Strassler solution, one expects the physics 
at the tip to be dominated by the anti-D3 branes sources, and in addition to have an extra D3- 
charge dissolved in flux, as expected from the brane-flux annihilation process [5]. As a result, we 
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need to consider how branes and fluxes contribute to the Maxwell charge. Recall that the Maxwell 
charge for D3-branes is defined as 

Qm x = Tj^-J ^\ (87) 



(4vr 



7' 11 



where we have indicated that the integral is over the T 11 slice at infinity. For the zero-th order KS 
solution, the Maxwell charge is not constant, it measures the running of the rank of dual gauge 
groups and thus is a rough measure of the degrees of freedom. However we are concerned only 
with the additional Maxwell charge which we will add to the background and thus we will abuse 
terminology and refer to this additional charge as Q^ x - In other words we have normalized the 
background Maxwell D3 charge to zero. 

As is clear from the Bianchi identity for F^ 5 \ Q^ x gets contributions from branes and flux 
combined 

djrO) = #(3) A F (3) + 5(x - Xi) (88) 

i 

where the second term comes from explicit D3-brane sources. Integrating flHTl) by parts yields 

Qm x = 9, + Of (89) 



where 

<lb 



(4tt 2 ) 2 
1 

2T2 



7111 

1 o 



(4tt 



F (5) , (90) 

H m A F ® (91) 



are the charges from the explicit D3-brane sources and flux respectively. Note that Q^ x is gauge 
invariant (unlike the Page charge) but is not necessarily quantized"]. 



6 Constructing the Anti-D3 Brane Solution 

As advertized in Section [TTT1 we have two strategies to look for the "anti-D3 brane in KS" solution. 
The first is to impose UV boundary conditions, perturb the KS solution, and try to match the 
expected boundary conditions in the infrared. The second is to impose the IR boundary conditions 
(which make the near-tip solution to be that of anti-BPS 3-branes in anti-Klebanov-Strassler), and 
to try to obtain the UV KS boundary conditions as a non-normalizable perturbation of the anti-KS 
solution. 

6.1 UV — > IR, Xi is non zero 

The UV boundary conditions we wish to enforce on our prospective new solutions are that Q^ x = 
N and that there should be a non- zero force on a probe D3-brane. In our ansatz the first boundary 
condition requires a careful analysis of 

jr 5 = (fcF + /(2P-F))vol T u (92) 
11 A very nice description of the various charges in type II string theory can be found in |29 | I30 | I5T]. 
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in the UV. In addition the brane probe calculation of Section |3] demonstrates that the force on 
a probe D3 brane depends only on X^ v ■ Since X[ R is proportional to X^ v (see Eq. (1701) ). 
demanding a non-zero force in the UV will imply 

X^ v ^ X{ R ^ . (93) 

In addition, we also demand very divergent UV non-normalizable modes to be absent, and hence 
set 

Y 5 UV = Y" v = 0. (94) 

These are the only UV conditions that we will impose so in principle all the other non-normalizable 
modes could be turned on. 

Now we need to extract from this abstract discussion the explicit UV and IR behavior for the 
modes in our deformation spaced. The total Maxwell charge is evaluated in the UV by demanding 

(06 - k)F + {ko - h)h + 2P05 = Qm x + 0(e" T/3 ). (95) 
Using (lT3"j) , f |94p and the UV expansions in section 14.3.31 we find that 



UV 



= 2PY 7 UV . (96) 

In the IR, we want the 0* to contain no divergent terms except those coming from q b . From 
f )90|) and ( TXQ|) we immediately find that, 

(0 6 - 5 )F O + (k - / o )0 7 + 2P0 5 ~ q b + 0(t) . (97) 

From the zero-th order expressions (1131) we see that 

F ~O(t 2 ), f ~O(T 3 ), k ~O(r) (98) 

so that disallowing divergences in (0 5 ,0 6 ,0 7 ) we find the constraint 

2P05 ~ q b + 0(t) . (99) 

Then on physical grounds we demand that in the IR the singularity in the Einstein tensor (due 
to the singularity in the warp factor) is commensurate with the singularity in the energy density 
from (recall that 04 is the warp factor): 

4> 4 ~\q b \/T + O(T ). (100) 

So this divergence is allowed since it is clearly coming from the anti-D3 brane sources. 

However we must eliminate other divergences. We can see from the equations for 2 and 3 
that this implies 

ylR = ylR _ q ( 1Q1 ) 

Note that unlike the integration constant X±, whose values in the UV and IR expansions are 
proportional to each other, all the other 13 integration constants differ the UV and IR by an 



2 Recall that in our notation (05, <pg, <fij) — (/i, fci, F\), the perturbations to (Jo, ko, Fq 
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additive constant, as reviewed in section fl~2l Hence, the fact that Y^ R must be set to zero implies 
that Y 2 UV will be a very nontrivial function of all the JQ, which is highly unlikely to be zero 
unless all the are zero and we have a BPS solution. Hence, any prospective anti-D3 solution 
will source modes that behave in the UV as 1/r 3 , contrary to what one might expect based on 
codimension-counting. 

The absence of divergences in 7 further implies that 

Yi* = ^{Xl* + Xl*) (102) 

and the absence of 1/ r-divergences in 6 implies 

= (103) 

ho 

Furthermore the absence of a log-divergent mode in 6 combined with the equations above also 
implies 

X{ R = (104) 
o 

The absence of a divergent term in 8 imposes a relation between X$ and the other variables: 

X IR = -j(5Xl R - X IR - 6Xl R ). (105) 

Furthermore, upon using ( 1103 j) one can see that the - term in 4 depends only on X[ R and Y 7 IR 
l 2 10 / 3 3 2 / 3 P_ rR 2 4 / 3 rp 4096P 4 



74 — - 

T 



( _ + ^ A ,, R)+y/fi + o(T) (io6> 



Upon using these conditions above, the expansions of 5 becomes 

5 = Y 7 IR + ^^(l) 2/3 Xl R + A(2/3Y/ 3 PXi R T 2 + +0(r 3 ). (107) 

Note however that the values for (Y 7 IR , X[ R ) have already been fixed in terms of N: firstly, ( |96i) 
fixes Y 7 UV in terms of the number of initial anti-D3 branes (N) put into the system, which is in 
turn related to Y 7 IR by an additive constant. Secondly, eqns. ffTOO]) and ffT06]) give X[ R in terms 
of Y 7 and q^. One can then use (I99T) and (11071) to obtain X( in terms of Y 7 IR . At this point we 
have computed all the IR quantities including q b in terms of iV and one might declare this to be 
the solution for anti-D3 branes in the KS background that we have been seeking. However, the 
singularity analysis is more subtle than that presented above. 

6.1.1 IR Singularities 

Recall that the NS and RR three forms in our ansatz are 



H3 = ^ 
F 3 = 



U k ~ f)9n A (9i A93 + 92A £4) + dr A {fgx A g 2 + k'g 3 A g A ), 
Fgt Ag 2 Ag 5 + (2P - P) g 3 A g 4 A g 5 + F'dr A (g x A g 3 + g 2 A g 4 ) , 
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and in the IR the unperturbed metric is regular and is given by 

ds 2 10 ~ axdsl + a 2 (^dr 2 + (gj + gj + ^g%) + ^{g\ + g%) J ■ 

for some numerical constants {c\,c>i). From this we observe that if 

(k - f) ~ a + 0(t) (108) 
then the energy density has a term that divergences as 

Hi ~ a l (109) 

Furthermore, if 

/ ~ 6 + hr + b 2 r 2 + 0(t 3 ) (110) 
then there is a singularity in the energy density of the form 

^-^ + ^- (HI) 

The singularities from cto and b 2 are somewhat more benign than that from 6 1; since their divergent 
energy density integrates to a finite action^ (recall that y/G ~ r 2 + 0(t a )) 

I J d 10 xVGHl\ < 00. (112) 

The 61 term in (11 101) however would result in a divergent action. Similarly, the IR expansion of 

4> 7 = co + Ci t + 0{t 2 ) (113) 
gives a singularity in the energy density of order 



and so we see that both the RR and NS three-form field strengths have potentially inadmissible 
singularities. 

Using the boundary conditions from section 16.11 we find that 

h-k = -^(i) 1/ V a + o(T) , (115) 

h = ^(D^r + OfA (116) 

The difference 05 — 06 has a nonzero constant in its IR expansion {clq 7^ 0), 05 has a nonzero 
quadratic term (b 2 7^ 0) and 07 has a nonzero linear term (ci 7^ 0), all of which are proportional 
to X-[ fi and lead to an energy density which diverges as t~ 2 . However, both the linear term in 5 



13 



We are very grateful to Igor Klebanov for sharing this observation with us. 



30 



and the constant term in (pj are zero (61 = cq = 0) and thus there are no singular energy densities 
that scale like r -4 , are hnce no infinite-action divergencies. 

Our analysis establishes that the only solution whose ultraviolet is consistent with that of 
would-be anti-D3 branes in the Klabanov-Strassler solution must have finite-action divergences in 
the infrared, both in the RR and the NS-NS three-form field strengths. The only way to eliminate 
these divergencies is by setting X[ R = 0, which in turn implies that X^ v = and is thus in 
contradiction with the UV boundary condition ( 1931) necessary in order to have a nonzero force on 
a probe D3 brane. 

Clearly, finite-action divergencies are better than infinite-action divergencies, but it is not clear 
whether in the context of the AdS-CFT duality even such finite-action divergencies should be 
allowed. The best example of a finite-action singularity that must be excluded from the spectrum 
of any consistent theory of two derivative gravity is the negative-mass Schwarzchild solution^ 4 ! 

Furthermore, this singularity does not appear a-priori to have any direct connection with the 
infrared physics of probe anti-D3 branes in the KS solution [5J, and the brane-flux annihilation 
that occurs via the polarization of these anti-D3 branes into NS5 branes wrapping an S 2 inside 
the S 3 : First, brane-flux annihilation results in the non-trivial division of Q^ x into qf) and 
this in turn is related to the IR and UV values of /, and not to the IR value of k — f. As such, 
this singularity could not have been predicted as a consequence of brane-flux annihilation. Second, 
the NS5 brane that mediates brane-flux annihilation couples magnetically with an NS-NS three- 
form field strength, while one of the fields whose energy density diverges in our solution is a RR 
three-form field strength (ci 7^ 0). Third, in the probe analysis of brane-charge annihilation, the 
amount of charge dissolved in fluxes depends nonlinearly on the amount of anti-D3 branes, while 
in our solution both the anti-D3 charge and the coefficient of the divergent terms are proportional 
toX( R . 

Hence, given the absence of a microscopic explanation for this singularity, and given the ab- 
sence of any criterion under which this singularity would be acceptable while the negative-mass 
Schwarzschild singularity would be not, it is fair to say that if this singularity had appeared all by 
itself in a fully-backreacted solution, it should have been deemed unacceptable. 

However, as we have seen in the previous section, this singularity does not appear all by itself: 
the infrared geometry also contains metric and five-form fields whose energy density computed in 
first-order perturbation theory diverges as r~ 4 . These fields have an obvious microscopic inter- 
pretation: they are sourced by a smeared distribution of anti-D3 branes, the metric and five-form 
divergences are commensurate, as typical for D-branes, and we expect this singularity to make per- 
fect sense in string theory. Hence, it is possible that upon taking into account the string-theoretic 
resolution of the singularity sourced by the smeared D3-branes, the the weaker finite-action sin- 
gularity may become benign. 



6.2 IR UV 

We now start from a solution that eliminates the IR singularities of the previous section, and argue 
that the UV expansion of this solution cannot match the UV expansion of the KS solution with 
none but the charge non-normalizable mode turned on. This problem is equivalent to starting 
from the BPS D3 branes in the IR and arguing that one cannot match the anti-KS solution in the 

14 We are very grateful to Gary Horowitz for pointing this out. 
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uv. 

If we are to match the boundary conditions near the anti-branes, Eqs. f ll0ip -( fT05l) . the fields 
5 and 6 must be equal, and using fjlOTj) this implies 

x[ R = X? v = 0. 

Note that in the absence of X± we have already solved exactly the equations (JTJ for the in 
Section Ei~Tl and thus the values of the Xi in the UV and in the IR must be the same. The infrared 
regularity equations in the previous section imply: 

x 5 + X 6 = X 4 = , 
Y 2 IR = Yj R = Yq R = , (117) 

and therefore 

X 8 = -P(X 5 - X r ) . 

Furthermore, since we do not want divergent and non-normalizable 1/r modes in the UV, 

X% = Xi = Xj = . 

However, in deriving the infrared boundary conditions in the previous subsection we have used the 
zero energy condition (IBlj) to express X 2 as a function of X 3 , X 4 , X 5 and X 7 , and setting it to zero 
combined with the other equations implies 

X 5 = 

and therefore all X's must vanish. Hence, if one tries to match all the correct IR boundary 
conditions and avoid un-physical non-normalizable modes, one must set all the to zero. 

One can now proceed exactly as in Section 14.31 to establish that <p 8 is a constant, and use 
equations ( |75l) . ( J76l) and ( I771) together with ( 111 T[) to show that 2 = 3 = 0, and that 1 must be a 
constant. The next step is to solve the system (1781) exactly, and match the infrared and ultraviolet 
integration constants. One finds that 



5 = ^sech 2 ^ ( cosh r(2V 5 r - V 5 sinh t + Y 7 -PY 8 t) 

+ sinh t(PY 8 - 2V 5 ) + V 5 r + V 6 + Y^ 

6 = - ^csch 2 1( - cosh r (2 Y 5 r + Y 5 sinh r + V 7 - PV 8 r) (118) 

+ sinh r(PV 8 - 2V 5 ) + Y 5 r + V 6 + V 7 ) , 

7 = — Y 5 coshr + (Y 5 t — V 6 )cschr . 
Comparing to the expansions in Sections 14.3.21 and 14.3.31 we get 

Y 5 = Y 5 UV = Y 5 IR + ^, Yq = Yq R = Yq V , Y 7 = Y 7 IR = Y 7 UV . (119) 

Requiring no divergencies in the UV and IR sets 

Y 5 = Y 6 = 0, (120) 
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Importantly, we see that the sign of Y7, which is proportional to the effective D3-brane charge 
(recall that Y% just corresponds to a shift in the dilaton), cannot change from the IR to the UV. 
This implies that this solution has positive D3 brane charge throughout Hence, the solution that 
preserves D3 boundary conditions in the infrared cannot match a solution with anti-KS boundary 
conditions in the UV and viceversa. 



7 Discussion 

With a view towards identifying the modes sourced by anti-D3 branes, we have computed the 
spectrum of linearized fluctuations around the Klebanov-Strassler solution which respect the sym- 
metries of the original background. Finding such solutions, with anti-D3 brane boundary condi- 
tions in the infrared and Klebanov-Strassler boundary conditions in the ultraviolet has been a key 
unresolved step in many constructions of de Sitter vacua in string theory. Much work has been 
done using this mechanism to show that large numbers of such vacua exist. 

We have found in our spectrum a single solution where some of the anti-D3 charge is dissolved 
in flux, supersymmetry is broken (as suggested in [5]) and there is a force on a probe D3 brane of 
the same form as the one computed in [13]. While this is all somewhat encouraging, the solution we 
find must necessarily have an infrared singularity, coming from the fact that the energy densities 
of both the RR and NS three form fluxes diverge quadratically at the origin. 

We have discussed this singularity in great detail in Section 16.11 and have seen that as far as 
first-order perturbation theory is concerned, this singularity is better behaved than other possible 
IR singularities - in particular its divergent energy density integrates to a finite action. However, 
even if the finiteness of the action makes this singularity more acceptable than others, it cannot 
be used as a criterion to establish that this singularity is physical. For example, the negative-mass 
Schwarzschild solution also has a finite-action singularity, and yet this solution is unphysical, and 
must be eliminated from the bulk spectrum in order for the gauge/gravity duality to make sense. 
Indeed, the breakdown of general relativity at a singularity does not come from the diverging 
action but rather from the fact that the derivative expansion in the effective action for gravity 
does not provide a good perturbation theory. Hence, divergent energy densities in the RR and 
NS-NS fields are as unphysical as a divergent Ricci scalar, and indicate that perturbation theory 
has broken down. 

We have also tried to find a microscopic interpretation for this singularity, to see whether the 
divergent fields can be explained as coming from D-brane or NS5-brane sources that one may 
expect to find in the infrared, and hence whether stringy corrections could resolve this singularity. 
We have argued that the existence of this singularity is neither a direct consequence of the expected 
brane-flux annihilation, nor of the polarization of anti-D3 branes into NS5 branes in the infrared. 
From such arguments, the most direct conclusion would be that this finite-energy singularity should 
be treated as other such singularities in gravity and must be removed. 

On the other hand there is one reason that we should persist with this solution and perform 
further tests on it despite the singularity. It may be that the finite-action singularity is not a 
signal of the breakdown of general relativity but rather a consequence. It is clear that one can 
never completely capture the physics of the solution near the anti-D3 branes by doing perturbation 
theory around the original KS background, since in that region the energy in the fields sourced by 
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the explicit anti-D3 brane sources diverges. This is similar to the fact that one cannot treat an 
electron in a background electric field as a perturbation, since near the electron its electric field 
always dominates. The energy density of the fields sourced by the explicit anti-D3 brane sources 
in the IR diverges quartically in the IR, whereas the energy density of the three forms diverges 
quadratically. It may be that the stringy resolution of the anti-D3 brane singularities, or even 
their full supergravity backreaction, will also resolve the milder finite-action singularity we have 
found and render the whole solution physical. 

It is clearly crucial to establish whether this finite-action singularity is indeed resolved this way 
and could be acceptable from the point of view of string theory, or whether it is not acceptable 
and must be removed. 

If the finite-action singularity is physical, than our analysis has found a first-order backre- 
acted solution that has anti-D3 brane charge at the bottom of the Klebanov-Strassler solution. 
The objects with anti-D3 brane charge source modes that in the ultraviolet can be treated as 
a perturbation of the Klebanov-Strassler solution, and it is very likely that this solution repre- 
sents the backreaction of the smeared system of anti-D3 branes considered by Kachru Pearson 
and Verlinde in [5]. Our analysis does not fully establish that these anti-D3 branes source only 
normalizable modes - to prove this one needs to relate the UV and IR integration constants, which 
is quite nontrivial |32j. If this is so, then this solution would represent the supergravity dual of 
the conjectured metastable vacuum of the Klebanov-Strassler field theory, and would be the first 
supergravity solution dual to a metastable vacuum of a four- dimensional theory. By exploring the 
vacuum expectation values of various operators, this solution would allow one to characterize very 
precisely the physics of this metastable vacuum, and would quite likely increase our understanding 
of gauge theory metastable vacua in general. 

On the other hand, if this singularity is not physical, our analysis establishes that the solution 
corresponding to backreacting anti-D3 branes at the bottom of the Klebanov-Strassler background 
cannot be treated in the UV as a perturbation around the BPS vacuum. Hence, these anti-D3 
branes source non-normalizable modes (other than the charge) whose energy is at least as strong 
as that of the original KS solution. 

While this conclusion might seem at first glance unexpected, it is is fact not once one realizes 
that some of the modes may have two-dimensional dynamics, and a perturbation is such modes 
does not decay at infinity. Indeed, as we have argued in the Introduction, all supergravity solutions 
dual to four-dimensional theories with a running coupling constant must have bulk modes that 
depend logarithmically on the radial direction, coming from the fact that the coupling constants of 
these theories run logarithmically with the energy. Hence, it is quite natural to expect the anti-D3 
branes to couple to those modes, and hence to source perturbations that grow logarithmically with 
the radius. It is not hard to see that the energy of such logarithmic modes is as strong as that of 
the unperturbed solution. 

Furthermore, this behavior would mimic exactly what happens when one tries constructing an 
MQCD metastable vacuum [I]. First, if one tries to obtain this vacuum by perturbation theory 
around the BPS one, one also finds modes that diverge in the infrared, and cannot match the 
correct boundary conditions. This is a first hint that the fully-backreacted metastable vacuum 
might not exist, and is confirmed by the subsequent backreacted calculation. Furthermore, if one 
tries to find the solution whose infrared has regular boundary conditions, one finds a solution that 
differs from the BPS one by a logarithmically-growing mode, and that cannot be obtained as a 
perturbation of the latter. This is again a reflection of the fact that this mode has two-dimensional 
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dynamics, and once excited in the infrared it does not decay at infinity. 

If one takes the analogy with the MQCD construction of metastable vacua one step further, 
one can also identify an obvious candidate for the full solution consistent with anti-D3 branes 
in the infrared and no finite-action singularity: the anti-D3 branes in the anti-Klebanov-Strassler 
solution. Much like the MQCD IR-compatible solution^, this solution differs from the BPS one by 
log-growing modes, it is not a perturbation of the latter and is BPS by itself. If this is indeed the 
solution for the back-reacted anti-D3 branes, this would imply that the sign of the charge dissolved 
in flux in the KS background must always have the same sign of the D3 branes at the bottom, and 
changes sign as one changes the orientation of these branes. 

If the finite-action singularity is excluded, or if it is allowed but the anti-D3 branes source modes 
that are nonnormalizable, this would strongly suggest that the dual gauge theory does not have a 
metastable vacuum. Though unexpected, this result would be consistent with the fact that a careful 
analysis of the vacuum structure of this theory [33] has not found such a vacuum. Furthermore, 
as we have explained in the Introduction, this would probably affect the KKLT construction of 
deSitter vacua in string theory, and the existence of a landscape of such vacua in general, which 
in turn may invalidate many string phenomenology and string cosmology constructions. 

It is fair to say that neither our analysis, nor any in the literature definitively show that the 
finite-action singularity is unphysical, and hence rule out the existence of KS anti-D3 metastable 
vacua. On the other hand, we have not been able to find any microscopic arguments or singularity- 
resolution arguments of why this singular behavior in the fields may be acceptable, which would 
establish that the solution we have constructed is the supergravity dual of the KS metastable 
vacuum, and implicitly prove that such a vacuum exists. 

Our analysis crystalizes and reduces the complicated problem of the existence of KS metastable 
vacua to the simpler problem of whether the finite-action singularity is physically acceptable or 
not. Establishing whether this singularity is physical would ideally come from a microscopic un- 
derstanding of the infrared physics or from a near-brane analysis of the fully backreacted, localized 
solution (for example extending the analysis of [SI])- Any developments on these issues would 
clearly be very exciting. 
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A Conventions 

We use the somewhat unsightly notation £ a [20] to indicate that we refer to £ a as defined in reference 
|20j . This is to try to keep straight the different conventions which appear in the various papers 
on the subject . The relation to the variables used in other papers is the following. The angular 

15 Given by Eq. (3.22) in @]. 
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variables e\, 62, €3, e\, e<i of [16] are related to Oj by 

ei = 75(^2 + #4) , e 2 = -^(-9i + 93) , ei = ^(# 2 + #4) , e 2 = ^(#1 + #3) , £3 = 9s ■ (121) 

Our radial variable r is related to the radial variables in [21 EH [201 IIZ1 EE2] by 

dr = dr^ = e 4pi rfrrj7] = e 4pi dr^ = e A{po+Pl) du^ = V6e 3po+Xo+4pi d\nr^ . (122) 

To zeroth order in perturbation this last relation reads in the UV r = e T//3 . 

Our metric functions x,y,p,A are related to q,y,p,Y used in [20] and a, b (y = and in the 
solution c = 0) in [12] (relations to the latter are only valid in the UV) by 

x = 3q + 2p[20] = -2a - ln6 , y = ygjg , 

p = Pgjg - 9 = ^(2a - & + ln(3\/6)) , A = Y = a + lnr . (123) 
The functions /, fc, F in the 3-form fluxes are related by 

/ = ~ 2p fm = 2 -Am = 6/c n2] ' 

k = -2Pk m = 2k m = 6k m , (124) 
F = 2PF M =2F m . 



36 



B Behavior of £ 

B.l IR behavior of £ 

The IR behavior of £ a is the following 



16 X^(§) 1/3 P 2 r 2 



6 = + ^ + G(r 4 ) 

8X 1 ' li (2)"'p ,„ 16 X' R (|) 1/3 P 2 t 

i» - ^ + ^ + (^&-m, t+0( ,, (125 , 

8X^(j) 1/3 P 2 X|* + Xf ^ ! 
6 = 1 ho + 2— 

. 4X 1 /i? (f) 1/3 P 2 X^+X^X X IR -AX IR ,^ 2N 
+ , _^_kl + j + _ r r + 0(r2) 

Is = X 1 * + \p(hX{ R - Xi R - 6X«) 

+p ^ 32X^(|) 1/3 P 2 + X^ + X 6 ^ \ ^ + 2 p(x „ _ ^3 + 0(r4) 



2 



V 45/i 



6 = "3 Xl + 3^ + 3 Xs T +0(T) 

.X[ R 16X^(1) 1/3 P 2 X 4 ^-P(X 5 ^ + X 6 ^ )- 
J + 3h + 3 

+ ( A X i R + 2X 2 /R - PX IR - ^Xi R + 3PX 7 ^ r 

,X[* 56X^(1) 1/3 P 2 X|«+jPjXj* + Xj*) ^ 
"18 15 ho 6 ! 



B.2 UV behavior of £ 

The UV behavior of £ a is the following 
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i } = Xf v 'e- 4 ^ 3 (3-12r) + O(e- 10T / 3 ) 

| 4 = y - Xf y (3 - 12r) e" 4 ^ 3 + O(e- 10r / 3 ) 

| 5 = X^ - Ix^e" 4 ^ 3 + O(e- 10T / 3 ) 

fa = \x^e r + {2t(X™ - X™) + X™ + X™ + \X™) e" T + 0^) 

f 7 = -^XfV+(2r(X^-XH-^ + ^ + |^7 aV )e- r + 0(e- 3T ) (126) 

+2p( - 2(Xf y - XfV + (3X 5 ^ - X^ - 4XDr + X 6 + X^e^ + 
| 3 = x 3 ^(e 2r -4r-e- 2r ) + — X 1 f/l/ e- 4T / 3 (-23 + 140r) + O(e- 10T / 3 ) 

6 = ((-^ + PXr)r + Xr)e T 

+ ^_^X 3 ^+5PXr-4PX^)r-Xr + ixr + P(X 5 C7y -X 6 C7y -2X 7 c/y )) e^ + G(e^) 
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